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Spectrometers 

This Invention relates to spectrometers and, 1n particular, 
to methods of multiplexing the signals received by 
5 spectrometers. 

In general, laboratory spectrometers receive signals from 
one radiation source and one sample, and display the results as 
a single spectrum of the sample. 

In many applications for sensing and process control 1t 1s 
10 desirable to use a single costly spectrometer to analyse the 
signals from many sources and/or samples. The simplest way to 
do this Is to present each signal sequentially to the 
spectrometer, and analyse the results separately. 
Alternatively, a multiplexing technique 1n which several signals 

15 are presented and analysed at once may be used. Multiplexing 
has two advantages, the total signal arriving at the detector Is 
greater than if the signals are observed sequentially and, for 
systems In which the noise Is Independent of the signal, the 
overall signal-to-noise ratio Is increased relative to that 

20 which can be obtained by sequential sampling. The disadvantage 
is that 1t Is necessary to determine the contribution of each 
source and sample to the total observed (composite) signal. 
This may be done by encoding the information carried by each 
signal 1n some appropriate way before recording the composite 

25 signal, and then using a suitable decoding process to recover 
the Individual signals. 

A variety of encoding and decoding schemes, which differ in 
their degree of difficulty of Implementation, are possible, 
Typically the observations Y 1 which measure the total signal 

30 presented to the decoder (data-handling system) can be related 
to the unknown individual signals Xj by a system of linear 
simultaneous equations which can be written in matrix form: 
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The unknowns Xj can be recovered by solving these equations. 
The equation (1) presents a unifying framework for 
discussing the design of various kinds of spectrometer. The 
simplest case 1s that of a dispersive Instrument, such as a 
conventional grating Infra-red spectrometer or continuous wave 
15 NMR. The values Xj are radiation intensities at various 
wavelengths, distinguished by a wavelength sorting device such 
as a diffraction grating. These intensities are recorded by the 
Instrument and output 1n the form of a spectrum , 
1»l,2,... t N. Thus In this trivial case, the coefficients Ajj 
20 comprise a unit matrix. 

More complicated, and closely related, examples are provided 
by Fourier transform Infra-red (FT-IR) and pulsed Fourier 
transform nuclear magnetic resonance (FT-NMR) spectrometers. In 
the former case, wavelength information (i.e. radiation 
25 intensities Xj at various wavelengths) is encoded by a Michel son 
Interferometer into a series of Intensities Y i sorted as a 
function of the position of a movable mirror 1n one arm of the 
interferometer. In the NMR case the encoding is achieved, not 
by some physical device, but by the precession and progressive 
30 dephasing of the nuclear spins within the sample. Intensity 
values Xj as a function of frequency are encoded and recorded as 
a series of induced voltages measured at successive time 
Intervals after the excitation of the sample. 

In both these cases, the matrix A is a unitary matrix of 
35 Fourier coefficients. Equation <1> can be solved by carrying 
out an inverse Fourier transform of the data Yj. 

An infinite number of other possible choices for the 
coefficients A^j can be envisaged. However, there are four 
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Important limitations on the possible choices: 

1. It must be possible to solve equation (1). The basic 
requirement 1s that the determinant of A should be 
non-zero. In order to obtain accurate numerical solutions, 

5 It 1s also desirable that A should not be ill-conditioned. 

2. It must be possible to realise the coefficients A^ 
physically in some way. The encoding process may be 
Implemented explicitly In hardware (e.g. the Hlchelson 
Interferometer), or It may be Implicit in the physics of the 

10 system (e.g. nuclear precession). The hardware 
implementation is typically easier 1f the array Is cyclic 
(I.e. each row Is obtained by shifting the previous row one 
element to the left or right). 

3. It Is desirable that the experimental errors or 
15 signal-to-noise ratios of the decoded Xj values should be at 

least as good as, or better than, the errors obtained when 
measuring the Xj values directly by some other technique. 
Further, the Improvement In s1gnal-to-no1se ratio should 
preferably be the same for all the Xy This Imposes a severe 
20 restriction on A. 

4. It Is desirable (though not essential) that the matrix A 
can be constructed and Implemented for any order N. Changing 
N should not require radical redesign of the instrument. 

We now consider points 2 and 3 above 1n more detail. 

25 Consider first a simple series of N measurements Xj of some 
quantity, such as the Individual weights of a group of objects. 
Suppose that the errors In these measurements have standard 
deviation a, and that the weights of arbitrary collections of 
these objects can also be measured, also with error a. It is 

30 possible to determine the weights of the individual objects more 
accurately than o by weighing collections of the objects in 
groups, rather than individually. The choice of objects in each 
group is made according to a suitable weighing scheme. Given N 
objects, N groups of these objects must be weighed according to 

35 a weighing scheme which may be represented by an N x N matrix of 
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zeroes and ones. The presence of a one In the (1,j)th element 
indicates the inclusion of object j in the 1th group, while a 
zero Indicates Its absence. In equation (1), the matrix A Is 
the weighing scheme and the quantities Y ? are the weights of the 
5 groups. For example, a suitable weighing scheme for 3 objects 
might look as follows:- 

110 1 

A» 101 (2 ) 
10 L Q 1 1 J 

A matrix of zeroes and ones 1s appropriate to a type of 
measurement 1n which the Individual measurements can be co-added 
into the group measurements, or omitted. In the weighing 

15 example, this is the situation when a spring balance Is used: 
each object can be put on the pan, or not, as required. 

If a beam balance Is used, we have three options for each 
object: use the left-hand pan, the right-hand pan, or neither. 
The weighing scheme for this experiment requires a matrix of 

20 zeroes, ones and minus ones. This clearly Illustrates that the 
nature of the experiment determines the possible values of the 
coefficients A^ which can be realised. 

Considerable theoretical work has been carried out to find 
suitable matrices of the above types, which satisfy the above 

25 limitations. The conditions are very restrictive, and suitable 
matrices known only for certain values of N. Confining our 
attention to matrices of zeroes and ones, the best solutions are 
the so-called Hadamard simplex matrices, or S-matrlces. The 
process of solving equation <1> Is then an Inverse Hadamard 

30 transform. Exactly one S-matrix (plus Its cyclic permutations) 
exists for each N of the form N = 4n-l , n = 1,2,3,... (Note that 
requirement 4 above is not satisfied. How we proceed when N Is 
not of this form will be dealt with below). Constructions for 
S-matrices of various orders satisfy requirement 3; the 

35 signal-to-noise ratio of each unknown Xj (computed by Inverse 
Hadamard transform) Is improved by an enhancement factor 



CN + 1) /N 
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2/N 2 

D 

relative to the s1gnal-to-no1se ratio of Xj measured directly. 

These principles have been used In the construction of 
Hadamard transform spectrometers (HTS). For example, Decker and 
Harwit (J.A. Decker, Or., and M. Harwlt, Appl. Opt., £. (1969), 

10 p. 2552) constructed a multiplexing dispersive Instrument using a 
multlsllt encoding mask at the exit focal plane of a 
conventional dispersive monochromator . This mask allowed more 
than one spectral resolution element (I.e. radiation at 
different wavelengths) to Impinge upon the detector at once; 

15 hence, the encoding mask multiplexed the dispersed radiation. 
The pattern of slits In the mask modelled the pattern of zeroes 
and ones In one row of an S-matrix. The mask was physically 
moved to allow a different combination of spectral resolution 
elements to impinge on the detector before each detector reading 

20 was performed. 

Previous applications of Hadamard transforms to spectroscopy 
have concentrated on obtaining multiplex and throughput 
advantages from spectrometers designed to look at one sample at 
a time. Me have now developed a Hadamard encoding scheme for 

25 multi-sample spectrometers which can examine several samples at 
once, a capability which finds particular application in such 
applications as on-line process control and monitoring. 

According to the present invention there is provided a 
spectrometer for the simultaneous examination of a plurality of 

30 samples comprising a source of radiation, primary encoding means 
to encode radiation from said source, sample receiving means to 
position said plurality of samples 1n the path of radiation from 
said source, secondary encoding means selectively to modulate 
the intensity of radiation directed to Individual ones of said 

35 plurality of samples by a predetermined amount and detector 
means to detect radiation received from said plurality of 
samples. 
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Referring to equation (1). The quantities Y ? , Xj are now no 
longer to be regarded as simple numbers, but Instead become 
vectors, Vj, Xj In a space with an arbitrary number of 
dimensions P. The P components of each of these vectors are to 
5 be regarded as radiation Intensities (or voltages, or any other 
kind of signal amplitude) measured at a series of P frequencies 
or times. Thus each vector represents a complete spectrum (or 
Interferogram, free induction decay, or any other Mnd of signal 
trace). For purposes of illustration, the vectors will be 

10 regarded as FT-IR interferograms 1n what follows. 

Composite interferograms t i are measured experimentally as 
additive combinations of some or all of the interferograms Xj 
contributed by individual samples. Thus the coefficients 
are either zero or one 1n this simplest version of the 

15 technique. The co-adding of the Xj 1s a process of vector 
addition, e.g. 

N 

?1 "I A n X* (4) 

in which the corresponding data 'in each component of the 
interferograms X-j are added to give one component of ? 1 . Each 
row of equation (1) represents a complete experiment in which 

25 signals from several samples are co-added to give a composite 
interferogram. Multiple scans may be performed in each 
experiment in order to improve signal-to-noise ratio by 
co-adding the scans (signal averaging). 

The elevation of equation (1) to a vector equation makes no 

30 difference to the previous arguments in favour of S-matr1ces as 
the best choice for the sampling matrix A. Provided the noise 
in the instrument arises only from the signal detection system, 
and is independent of the numbers of sources and samples 
contributing to the composite signal, the use of S-matrices as 

35 the signal encoding scheme will lead to an improvement in 
signal-to-noise ratio (equation (3), with N the number of 
samples). 
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In the above example where the signals are regarded as FT-IR 
Interferograms, the experimental data (the N interferograms 
are doubly encoded signals. To recover the spectra of the 
individual samples, we need to perform (1) an Inverse Hadamard 
5 transform to recover the individual Interferograms, followed by 
(11) apodlsation and Inverse Fourier transforms to obtain the 
spectra. 

This technique exhibits advantages arising from multiplexing 
and signal averaging, and, in addition, It has the important 
10 advantage that only one expensive spectrometer or other signal 
detection system 1s required to measure multiple signals. If 
the radiation sources are Inexpensive, the performance can be 
Improved by employing one source per sample, thus Increasing the 
total signal arriving at the detector. 
15 The Invention will be particularly described with reference 
to the accompanying drawings, in which: 

Figures 1-3 represent basic multiplexing design 

schemes for different kinds of 
spectroscopic and other signal analysis 
20 Instrumentation. 

Figures 4a-d show four spectra composed of one or more 

Lorentzian lineshapes with Gaussian noise 
superimposed. 

Figures 5a-g show composite shapes obtained by 

25 combining the four individual spectra 

(the Lorentzian lines, not the noise) and 
then adding noise with the same standard 
deviation as in Figure 4; and 
Figures 6a-d show demultiplexed spectra calculated In 

30 accordance with one mebodlment of the 

invention. 

Some basic designs for multiplexing spectrometers are shown 
in Figures 1 to 3. All the schemes consist of one or more 
sources S, a primary encoder PE , a secondary encoder SE, an 
35 arbitrary number N of samples Sj and a detector D which converts 
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the signal to a digital record for analysis. In Figures 1 to 3 
only one source 1s used. In optical systems this would be a 
lamp or laser but could equally be a microwave generator or 
other source. 

5 Referring now to Figure 1 of the drawings, the energy from a 
source S 1s passed through a primary encoder PE. The primary 
encoder acts as a wavelength sorting device. It may be a 
diffraction grating, a simple filter system or a Fourier 
transform device such, as a Michel son Interferometer. 

10 Alternatively, for example in Fourier transform dielectric 
spectroscopy or NMR, It may be a pulse generator. The general 
function of the primary encoder is to ensure that each 
wavelength arriving at the detector is encoded in some way, to 
make it distinguishable from all the other wavelengths that may 

15 arrive. 

In Figures 2 and 3 the primary encoder 1s placed after the 
sampling region R; the total signal is then wavelength encoded 
after the radiation passes through the samples. In optical 
spectroscopy this means that some signal, already attenuated by 

20 absorption by the samples, acts as a source for the spectrometer. 
In Figure 2 multiple sources are used. This ensures that a 
larger signal energy will impinge on the detector. This 
configuration also has the advantage of improving the robustness 
of the system since the failure of a single source will not 

25 result in total signal loss but only loss of signal from one 
sample. 

A feature of all three configurations is the multiple 
sampling channels which are co-added to present a single signal 
to the detector. In general for optical systems fibre optics or 
30 light guides will be needed for at least part of the system, 
since the samples may be physically separated. Such devices are 
used throughout the system and transduction to an electric 
signal made at the detector stage. 

The secondary encoder SE is the device by which the signals 
35 are encoded by the coefficients A^j in equation 1. It modulates 
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the Intensities of the signals In the separate sample channels 
as required. If S-matr1x coefficients are used (Ajj ~ 0 or 1), 
then the secondary encoder takes the form of a series of 
computer-controlled switches, which may be mechanical shutters 
5 or a bank of electro-optic switches such as an LC-OSA, for an 
optical system; or a series of electronic gates, If the signals 
are electrical In nature. Two alternative positions for the 
secondary encoder are shown in each Illustration. 

Figures 1-3 represent basic multiplexing design schemes for 

10 different kinds of spectroscopic and other signal analysis 
instrumentation. In each case, multiple signals from several 
samples are encoded using successive rows of an S-matr1x to 
regulate the opening and closing of the switches In the 
secondary encoder. One problem remains: It was stated 1n the 

15 previous section that an S-matrlx only exists when the matrix 
order N Is of the form 4n-l, for n=l,2,3... How should we 
proceed when the number of samples does not satisfy this 
condition? The solution will be Illustrated with reference to 
an experiment using four samples. The smallest suitable 

20 S-matrix has N=7. Equation (1) takes the form, 
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Since samples 5, 6 and 7 do not exist, the vectors X 5 , X 6 
and X 7 are null vectors, or arrays of zeroes in the computer 
sense. At the experimental stage, we only need to employ the 
coefficients in the 7 x 4 submatrix on the left-hand side of A, 
35 That Is, we perform seven separate experiments to measure seven 

composite signals ? A Y 7 using the seven rows of the submatrix 

as the schedule for controlling the four secondary encoder 
switches or shutters. Thus, the seven composite signals are 
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composed from the four Individual signals as follows:- 



10 
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20 Since we are only Interested 1n 5^ X 4 we only need to 

employ the coefficients In the 4 x 7 submatrix at the top of the 

Inverse matrix A" x In the calculation. 

An important question Is whether the performance of seven 

separate experiments on just four samples In this way 1s more 
25 efficient than doing four separate experiments on the Individual 

samples. The answer, surprisingly, is yes, as will now be 

demonstrated. 

Suppose first that we perform one scan to acquire each of 
the seven composite signals. The slgnal-to-noise enhancement 

30 (the multiplex advantage) Is given by equation (3) with N=7, 
i.e. 1.512. Alternatively, we might perform four experiments on 
the individual samples with two scans per sample. The 
signal-to-noise enhancement obtained from signal averaging of 
two scans is /2, or 1.414. Comparing an enhancement of 1.512 1n 

35 seven scans with one of 1.414 in eight scans, we see a small but 
definite advantage In favour of the multiplexing scheme. 

Table 1 lists the expected signal-to-noise enhancements with 
sequential multiplexing sampling schemes, for an instrument with 
up to 15 sample channels. We see that for 3 or more samples, 

40 there is always an advantage in using a multiplexing scheme. In 
each case, the multiplex advantage is more than sufficient to 
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offset any redundancy Implied by the performance of more 
experiments than there are samples. 

A numerical simulation suffices to demonstrate that the 
above Ideas work 1n practice. Figures 4a-d show four spectra 
5 composed of one or more Lorentzlan Uneshapes with Gaussian 
noise superimposed. The exact nature of these spectra 1s 
Irrelevant to the discussion, but It should be noted that the 
standard deviation of the noise In each data set was the same. 
The spectra are plotted with a normalised vertical scale so that 

10 the most Intense peaks In each spectrum have the same height. 
Figures 5a~g show the seven composite shapes obtained by 
combining the four Individual spectra (the Lorentzlan lines, not 
the noise) according to equation (6), and then adding noise with 
the same standard deviation as In Figure 4. Finally, Figures 

15 6a-d show the demultiplexed spectra calculated according to 
equation (7). Comparing these with Figures 4a-d, an improvement 
(albeit small) in slgnal-to-noise ratio is evident 1n each 
case. The essential condition for obtaining this improvement is 
that the noise In the composite spectra 1s independent of the 

20 number of contributions. 

The Inverse Hadamard transform represented by equation (7) 
is a trivial mathematical process, requiring negligible computer 
time, because the number of samples (^the order of the matrix) 
will not generally be large. 



Number of 
channel s 



Order of 
S matrix 



Table 1 

Signal /noise 
enhancement 
for sequential 
sampling 



Signal/noise 
enhancement 
for S matrix 
sampling 



2 
3 



3 
3 



1.414 (4> 
1.0 (3) 



1.06 
1.06 



(3) 
(3) 



10 



4 

5 
6 
7 



7 
7 
7 
7 



1.414 (8) 
1.414 (10) 
1.414 (12) 
1.414 (14) 



1.512 (7) 

1.512 (7) 

1.512 (7) 

1.512 (7) 



15 



8 
9 
10 
11 



11 
11 
11 
11 



1.732 (24) 

1.732 (27) 

1.732 (30) 

1.732 (33) 



1.809 (11) 
1.809 (11) 
1.809 (11) 
1.809 (11) 



20 



12 
13 
14 
15 



15 
15 
15 
15 



2.0 (48) 

2.0 (52) 

2.0 (56) 

2.0 (60) 



2.066 (15) 
2.066 (15) 
2.066 (15) 
2.066 (15) 



25 The figures in brackets represent the total number of scans 
required to achieve the signal to noise enhnacement shown in the 
col umn 
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Claims 

1. A spectrometer for the simultaneous examination of a 
plurality of samples comprising a source of radiation, primary 
encoding means to encode radiation from said source, sample 

5 receiving means to position said plurality of samples In the 
path of radiation from said source, secondary encoding means 
selectively to modulate the Intensity of radiation directed to 
Individual ones of said plurality of samples by a predetermined 
amount and detector means to detect radiation received from said 
10 plurality of samples, 

2. A spectrometer for the simultaneous examination of a 
plurality of samples as claimed in claim 1 wherein a source of 
radiation is provided for each of said plurality of samples. 

3. A spectrometer for the simultaneous examination of a 
15 plurality of samples as claimed in claim 1 wherein said said 

primary encoding means is positioned between said sample 
receiving means and said detector means. 

4. A spectrometer for the simultaneous examination of a 
plurality of samples as claimed in claim 1 wherein said said 

20 primary encoding means acts as a wavelength sorting device. 

5. A spectrometer for the simultaneous examination of a 
plurality of samples as claimed in claim 4 wherein said said 
primary encoding means 1s a diffraction grating. 

6. A spectrometer for the simultaneous examination of a 
25 plurality of samples as claimed 1n claim 4 wherein said said 

primary encoding means is a Michelson interferometer. 

7. A spectrometer for the simultaneous examination of a 
plurality of samples as claimed in claim 4 wherein said said 
primary encoding means is a pulse generator. 

30 8. A spectrometer for the simultaneous examination of a 
plurality of samples as claimed in any one of the preceding 
claims including secondary encoder means comprises means to 
modulate the intensity of a plurality of signals Xj wherein, for 
a plurality of observations Yj such that 
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wherein said coefficients have the value 0 or 1. 



9. A spectrometer for the simultaneous examination of a 
plurality of samples as claimed in claim 8 wherein said 
secondary encoder comprises a Michel son Interferometer. 

10. A spectrometer for the simultaneous examination of a 
15 plurality of samples as claimed 1n claim 8 Including means for 

detecting precession and dephasing of nuclear spins within a 
sample. 

11. A spectrometer for the simultaneous examination of a 
plurality of samples as claimed 1n claim 8 wherein said 

20 secondary encoder includes a plurality of switches. 

12. A spectrometer for the simultaneous examination of a 
plurality of samples as claimed in claim 11 wherein said 
switches comprise a bank of mechanical shutters. 

13. A spectrometer for the simultaneous examination of a 
25 plurality of samples as claimed in claim 11 wherein said 

switches comprise a bank of electro-optic switches. 

14. A spectrometer for the simultaneous examination of a 
plurality of samples as claimed in claim 11 wherein said 
switches comprise a plurality of electronic gates. 
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